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ABSTRACT 

The ultraluminous X-ray source (ULX) Holmberg II X-l has been observed 
over 4 months in 2009/2010 by the Swift observatory. The source luminosity 
varied by a factor of up to 14, reaching a maximum 0.3-10 keV luminosity of 
~ 3.0 x 10 40 erg s _1 . The spectral properties do not vary much over these 4 
months, with only a slight monotonic increase of the hardness ratio with the 
count rate. This means that the erratic flaring activity of the source is not 
associated with spectral changes, as seen in other ULXs. Conversely, comparison 
with data obtained by Swift in 2006 shows a completely different picture: while 
at a luminosity also seen in the 2009/2010 data, the source appears with a hard 
spectrum. Thus, it appears that, as in Galactic black hole binaries, spectral 
states in this ULX are not determined only by the X-ray luminosity. 

Subject headings: accretion, accretion disks - black hole physics - X-rays: binaries 
- X-rays: individual (Holmberg II X-l) 

1. Introduction 

Ultraluminous X-ray sources (ULXs) are point-like X-ray sources located outside the 
nuclei of galaxies, with luminosities Lx > 3 x 10 39 erg s -1 . These apparent luminosities, 
assuming isotropic emission, are above the Eddington limit of a 20 M Q black hole. If the 
Eddington limit applies to these objects then the most likely e xplanation is accretion onto 



intermediate-mass black holes (IMBHs) of mass ~ 10 2 — 10 4 M (IColbert h Mushotzkvlll999 
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Makishima et al.ll2000l ). Alternatively, ULXs may represent a class of super-Eddington emit- 
ters and probe a regime in which the accretion r ate is much higher than that seen in Galactic 
black hole binaries (GBHBs; see iRobertsI 120071 for a recent review). The reality is probably 
more complex than a single answer, with some ULXs being possibly slightly more massive 
than stellar mass black holes (Mulx < 100 M Q ) but still experiencin g super-Eddington ac- 
cretion to some degree. In addition, mild beaming may be present (IPoutanen et al.l 120071 ; 
King|l2009h . 



R ecently, the Swift obse rvatory has been used to follow up some of the most luminous 
ULXs (IKaaret fc Fengl 120091 ) and has allowed studies that span a few months to a year with 
many observations, compared to the previous deeper but sparse Chandra/ XMM-Newton ob- 
servations. The results from these observations have shown that their flux varies significantly 
by factors of 5-10 on timescales of days to weeks ( Kaaret fc Fengl 2009 ). In NGC 5408 X- 
1, a 115 day orbital period has been suggested (jStrohmayerl 120091 ) . Somewhat unexpect- 



edly, most of the ULXs monitored with Swift to date (Holmberg I X X-l, NGC 5408 X 



1, and NGC 4395 X-2 ) do not show pronounced spectral changes (IKaaret &: Fengl 12009 



Vierdavanti et al.ll2010l) . The exception is the Hyperluminous X-ray source (Lx > 10 41 erg s 



Gao et al.l 120031 ) HLX-1 (IFarrell et al.l l2009h which shows X-ray variability (by a factor of 
> 20) but associated with spectral variability (IGodet et al.ll2009l ). 



Here, we report on a set of new Swift observations aimed at studying Holmberg II X- 
:eafte r) . This source is a luminous ULX (Lx up to 2 x 10 40 erg s _1 ; e.g. 
2009^. located in the dw arf galaxy Holmberg II, at an estimated distance 



Feng fc Kaaret 



d = 3.39 Mpc dKarachentsev et al 



20021). In optical wavelengths, the ULX is surrounded by 



the "foot nebula" (jPakull fc Mirionill2002l ) which is powered by reprocessing of X-rays and 
has enabled demonstration, using the He IIA4686 nebular recombination line, that the true 
X-ray luminosity is close to that inferred from the X- ray flux assuming isotropic emission. 



thus excluding strong beaming effects for this source (iPakull fc Mirioni 



2004 ). This has been confirmed by an infrared study ( Berghea et al. 201oh . using the [O IV] 



2002; Kaaret et al. 



25.89 /im emission line. This ULX has also been s 


tudiee 


. extensively in the X-ray band 


(Zezas et al.lll999; Mivaii et al. 2001; Dewanean et al. 


2004; 


Stobbart et al. 


2006: Goad et al. 


2006 


Feng & Kaaretll2009: Gladstone et al.ll2009: Kaiava & Poutanenl2009 


Caballero-Garcia & Fabian 


2010 


1 and has been seen with 0.3-10 keV luminosities ranging from a few 10 39 to ~ 2 x 



10 



10 



erg s" 



-3 - 



2. Observations and data analysis 

We obtained 68 observations under Swift program 90217 (PI: Kaaret) and also analyzed 
6 observations obtained under program 35475. All these observations were carried out using 
the X-ray Telescope (XRT) instrument in its photon-counting (PC) mode. The first data 
set is contiguous and spans 131 days from 2009 December to 2010 April and the second one 
comes from observations executed in 2006. We retrieved level 2 event files from all these 
observations. Thus, we used the default data screening parameters as described in the XRT 
user's guidfl Each observation is composed of one or more snapshots. We analyzed each 
snapshot separately, rejecting the ones with an exposure time below 100 s. We were left with 
124 snapshots. 

We extracted source counts from a circular region with a radius of 20 pixels (~ 47", 
corresponding to 90% of the point-spread function at 1.5 keV) and background counts from 
an annulus with an inner radius of 50 pixels (~ 118") and an outer radius of 120 pixels 
(~ 283"). The same regions were chosen to extract the spectra. The count rates and spectra 
were corrected for the loss of flux due to bad pixels and bad columns. For this, an exposure 
map was generated for each snapshot and used to create an auxiliary response file (ARF). 
The light curves were binned at 1 day averages to increase statistics. We also calculated a 
hardness ratio, defined as the ratio between the net count rate in the 1.5-10 keV band versus 
the 0.3-1.5 keV band. 

We co-added spectra using three different count rate ranges (see Table [1]) and keeping 
the 2006 data separate, using the addspec ftooiJE Thus, we ended up with four co-added 
spectra, which we refer to throughout the Letter as spectra 1, 2, 3, and 4 (Tabled]) where 
spectrum 2 is the spectrum for the 2006 data. The ARFs were co-added separately using 
addarf and weighted accordingly to their counts. Response matrix files (RMFs) were taken 
from the calibration database according to the date of the observations. Concretely, we 
used the RMF swxpc0tol2s0_20010101v011 . rmf for the 2006 observations and the RMF 
swxpc0tol2s6_20070901v011 .rmf for the 2009/2010 observations, grppha was used to bin 
the spe ctra with at l east 20 counts in each bin. Finally, we fitted the spectra using xspec 



12.6.0 flArnaudlll996h . 



*Can be found at http://heasarc.nasa.gov/docs/swift/cinalysis/ 

2 



2 We used FTOOLS 6.9, the latest version being available at |http :7 /heasarc . gsfc .nasa.gov/ftools/ 
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Results 



The X-ray light curve of HoII X-l (Figure [T]) shows significant variability, with count 
rates ranging from 0.02 to 0.28 counts s _1 . We calculated a periodogram from the 0.3-10 
keV light curv e from the 2009/20 1 data set, covering 131 days starting at MJD 55166, using 
the method of iHorne fc Baliunasl (119861 ) for periods in the range of 4-65 days. We found no 
significant peaks in the resulting power spectrum. 

The hardness/intensity diagram (HID) is shown in Figure|2j The data from the 2009/2010 
data set show a slight trend where the hardness ratio increases with the count rate. This is 
confirmed by a \ 2 test when we test for the null hypothesis of constant hardness. We get an 
average hardness of 0.34 ± 0.01 and x 2 — 150.6 for 54 degrees of freedom (DoF), ruling out 
constant hardness. A linear fit to the 2009/2010 data shows a good correlation between the 
hardness ratio and the rate, with a slope of 0.98 ± 0.10 and a corresponding \ 2 f° r the fit of 
52.2 for 53 DoF. 

However, the 2006 observations show a clear offset in hardness with a mean value of 
0.68 ± 0.03, at a mean count rate of 0.159 ± 0.004 counts s" 1 . At the same count rate, the 
mean hardness ratio is 0.38 ± 0.04 in the 2009/2010 data. Clearly, the hardness ratio is not 
solely a function of the count rate. 

We start the analysis of the spectra (Table [1] and Figure [3]) by using the most simple 
model, i.e., an ab sorbed power-law continu um. We fix the galactic H I column density to 



3.4 x 10 cm (iDickey fc Lockmanl |1990[ ) and allow for intrinsic absorption. The x °f 



the individual fits are statistically acceptable, except for spectrum 1. Spectra 1, 3, and 4 
show a consistent extragalactic absorption with a value of ~ 0.15 x 10 22 cm -2 and spectrum 
2 shows a lower absorption with 0.05 x 10 22 cm" 2 . The intrinsic colum n density is thus 



in reasonable agreement with previous studies (e.g.. iFeng fc Kaaretll2009l ) except the lower 
value for spectrum 2. Using the multicolor disk (MCD) blackbody model alone leads to a 
very bad fit, with a reduced \ 2 above 2.0 for all individual fits. 

Adding an MCD component to the power-law model improves the fits significantly in 
spectra 1 and 4, with A% 2 = 22.5 and 19.4 for 2 additional DoF. The inner-disk temperature 
is consistent within the errors between the different spect ra, with a value of ~ 0.20 keV, 
in agreement with results reported by other authors (e.g., IFeng fc Kaaretll2009f ). After we 



fitted the power-law index separately, we found that spectra 1 and 4 had the same index 
within the errors (Ti = 2. 80^22 •• -^4 = 2.57^Q'}g), consistent with the small evolution of 
the hardness ratio. Spectrum 3 which has fewer counts shows a compatible index but with 
a larger error (r 3 = 2.20^Q gg). Considering the relatively small variation in hardness ratio 
over the 2009/2010 data, which may be due to a changing fraction of flux in the soft thermal 
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component, we decided to tie these three indices together. The fit statistics are similar with 
X 2 /DoF = 380.2/375 when the photon indices are independent and x 2 /DoF = 384.4/377 
when not. We note that the "bump" we see in the residuals of the fit (Figure [3]) between 
2 and 5 keV is likely to be a real feature in the spectra (as seen in a recent XMM-Newton 
spectrum; J.J.E Kajava et al., in preparation). This "bump" is most apparent in spectrum 
1, where the reduced x 2 °f the fit is ~ 1.2. More complicated models should be able to 
model that feature but here we are limited by statistics. 

The power-law plus MCD model has been proven to be very effective in the study of black 
hole X-ray binaries and has been used extensively wit h the moderate signal-to-noise spectra 



nole A-ray binaries and nas been used extensively wrt n tne moderate signal-to-noise spectra 
of ULXs, leading to claims of the presence of IMBHs faaaret et alil2003l : iMiller et al.ll2003h 



due to the presence of a soft excess. But this excess may also arise from the emission of an out- 
flow at a large distance (s pherization radius) from the central source fjPoutanen et al.l 12007 
Kajava &: Poutanenll2009l ). Recently, some auth ors have shown that m ore physical models 



could be applied to high-quality spectra of ULXs. IStobbart et al.l ( 120061 ) and lGladstone et al. 



(120091 ) use a Comptonization model that reflects the coronal emission (and t hus replaces the 
power law) in addition to a disk model. ICaballero-Garcia &: Fabianl (120101 ) have applied a 
model based on a ionized disk reflection and power-law continuum. We a ttempted to ap 



ply on e of the Comptonization models (diskpn + comptt, in xspec) used by I Gladstone et al. 



(120091 ) to our spectra but the fits suffer from huge uncertainties in the fitted parameters. 
We also summed all the 2009/2010 spectra but were still unable to constrain the model 
parameters. 

It is apparent from the spectra (Figure [3]) and spectral parameters (Table [T|) that there 
are significant changes that are not only a function of the luminosity, as already seen in 
the HID. Spectra 1, 3, and 4 share the same power-law index (r ~ 2.6) and inner-disk 
temperature (T in ~ 0.20 keV) within the errors although there is a factor of ~ 5 difference 
in luminosity between spectra 1 and 4. The data from 2006 (spectrum 2) definitely show a 
harder power- law index (r ~ 1.8), but with the same T m and a luminosity close to spectrum 3. 



4. Discussion 



If we attempt a comparison with the spectral states of GBHBs, the power-law index of 
spectra 1, 3, and 4, which are typical of what has been seen during previous observations 
(r ~ 2.6) wo uld be consistent with some kin d of soft state. The thermal-dominant state 
as defined by iMcClintock fc Remillardl (120061 ) implies a disk-flux fraction above 75% and a 
power-density spectrum with no quasi-periodic oscillations (QPOs) or very weak features. 
The timing properties would be consistent here since the X-ray variability of HoII X-l is weak 
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Table 1: Spectral Fit Parameters 
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0.05 
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0.04 



Power-law + disk blackbody with Ti 
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+ 0.13 
-0.09 
+ 0.06 
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Note. — a Spectrum index used in the text ; b Count rate range used to create the spectrum ; c External 
absorption column ; d Power-law photon index ; e Inner-disk temperature ; f Absorbed flux (0.3-10 keV) ; 
s Unabsorbed luminosity (0.3-10 keV) for D = 3.39 Mpc ; h Fraction of the total unabsorbed flux (0.3-10 
keV) in the disk component ; 1 x 2 and degrees of freedom. All errors are at the 90% confidence level. The 
total number of counts in each of the four combined spectra is respectively 2314, 1940, 943, and 5368 counts 
for spectra 1, 2, 3, and 4. 



on short timescales f lGoad et al.ll2006l ). but the fraction of the flux in the disk component 
is below 50 % (with the caveat that the errors on the fractions are huge). Thus, it appears 
that in 2009/2010 the ULX would be classified as being in the steep power-law state where 
the power-law component contributes more than 50% without the presence of QPOs. 

Spectrum 2 shows a power- law index (r ~ 1.8) that would be indicative of the hard 
state. The disk component contributes ~ 14% of the total X-ray flux which would be 
consistent with the fraction seen in GBHB in this state. However, the temperature is the 
same, within errors, in all the spectra. This is unlike the behavior in GBH B, where the 



disk blackbody usually has a temperature well below that in the soft states ( iBellonil 12001 



Rodriguez et al.l 120031 : iRemillard Sz McClintockl 120061 1 . The properties of the soft thermal 
component in the spectra of HoII X-l are unlike the behavior expected for a disk blackbody. 
Thus, interpretation of the spectral states of HoII X-l in terms of the canonical states at 
GBHBs in which the thermal emission arises from an accretion disk is not valid. 

This particular ULX was only seen once in such a hard-like state (and at a similar 
luminosity than the Swift 2006 data; iMiyaji et al.ll200ll ) although it has been observed many 
times. Looking at the hardness ratios, it appears that HoII X-l stayed in this state during 
the three Swift observations executed in 2006 although they were separated by ~ 40 days. 
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Since HoII X-l has been caught only once before in such a state and the recent coverage 
of ~ 130 days does not show such a hard power-law index, interesting questions are raised 
such as does this "transi tion" happen rarely an d how long does it last. We note that the 
evolutionary sequence of iGladstone et al.l ( 120091 ) using Comptonized models is constructed 
as a function of luminosity and does not include the behavior seen here. Also, the constancy 
of the parameters of the disk component, discussed above, would require the corona to vary 
so as to cancel out the changes in the true disk temperature, which seems unlikely. It is 
also worth noting that the few persistent GBHBs behave quite differently compared to the 
transient ones. Indeed, objects like LMC X-l, LMC X-3, or Cyg X-l sp end most of their 
time in either the soft state or the hard state with only a few transitions (IBellonil |20 lOL and 
references therein). From a strictly spectral hardness point of view, HoII X-l also seems 
to spend most of its time in a soft state. Other ULXs like H olmberg; IX X-l seem to be, 
instead, locked in some kind of hard state (IKaaret fc Fengn2009i ). We also not e that the ULX 



NGC 1313 X-l was seen in two different "states" at a similar luminosity (IFeng fc Kaaret 



20061 ). hence showing an analogous behavior to HoII X-l. 



On the other hand, the 2009/10 light curve reveals that HoII X-l experiences rapid flares 
that are not tied to spectral changes (i.e., the hardness ratio does not vary much). These 
flares occur on timescales shorter than the spacing between the observations, < 2 days with 
a mean amplitude of ~ 0.15 counts s _1 (absorbed flux/unabsorbed luminosity amplitudes of 
~ 4.5 x 10 -12 erg s _1 cm _2 /1.2 x 10 40 erg s _1 ). If we consider all the snapshots without any 
additional binning, there are two cases that are indicative of a consecutive flux rise and decay 
in less than half a day each. Other ULXs also show large X-ray variability without significant 
changes in their spectral hardness (IKaaret fc Fengil2009l ). T he irregular flaring ac tivity that 
we see here seems very similar to that of Holmberg IX X-l (IKaaret fc Fengil2009l ) and could 
be also compared to the variability of the GBHB GRS 1915+105, which exhibi ts rapid flares 



on t imescales of minutes) that are associated with specific variability patterns ( iBelloni et al. 



200oh and interpret ed as an instability in the radiation pressure dominated part of the disk 



( IJaniuk et allboOOf ). 

We have shown that the ULX HoII X-l demonstrates X-ray variability by an average 
factor of ~ 5 on timescales of days to months. Over 4 months in 2009/2010, its spectral 
properties did not vary much, with only a slight increase of the hardness ratio associated 
with an increase in the count rate. During that interval, its spectrum was soft with a steep 
power-law component and a soft thermal component, similar to tha t seen during a number of 
previous observations (e.g.. lGoad et al.ll2006t IFeng fc Kaaretll2009l ). Only the fraction of flux 
in the soft thermal component has been seen to vary as a function of X-ray luminosity. This 
quasi- absence of spectral va riability over large changes in flux has been observed in other 
ULXs (IKaaret fc Fengn2009f ). But if we add the 2006 data, the HID becomes non-monotonic: 
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the hardness ratio is higher and correspondingly the spectrum shows a harder power law, 
consistent with the hard state, although its luminosity is close to two of the other spectra. 
The thermal component appears at a temperature close to that seen in the soft spectra, 
which argues strongly against interpreting that component as due to disk emission. Thus, 
while HoII X-l appears to exhibit hard versus soft states similar to GBHBs, interpretation of 
the thermal component in terms of GBHB states is likely invalid. Even though our spectral 
analysis is limited by the use of simple and somewhat empirical models, there are spectral 
changes that are clearly not asso ciated with flux changes. This behavior means that, as 
with GBHBs Jrloman et alJl200lh . at least two parameters determine the accretion state of 



this system and that the accretion rate is not the only variable responsible for spectral state 
changes. 

To understand the physics involved, deeper spectra at different count rates and hardness 
ratios are needed. The unprecedented coverage possible with Swift has led to new insights 
into the behavior of ULXs. Combined use of Swift with larger X-ray telescopes could lead 
to a better understanding of the evolution of the spectral properties of ULXs. 
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Fig. 1. — X-ray light curve in the 0.3-10 keV band for Holmberg II X-l (top) and corre- 
sponding hardness ratio (bottom). 
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Fig. 2. — Hardness/intensity diagram for Holmberg II X-l. Each individual observation from 
the 2009/2010 observing campaign is plotted as a green bar. The blue diamonds represent 
the 2006 data. The red dashed horizontal lines show the count rates ranges used to co-add 
the spectra. In each count rate range, we calculated the average hardness ratio, shown by a 
red cross for the 2009/2010 data or a red diamond for the 2006 data. 
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Fig. 3. — Swift XRT best spectral fit using an absorbed multicolor disk plus power-law 
model. From bottom to top: spectrum 1 (black), spectrum 2 (red), spectrum 3 (green), 
and spectrum 4 (blue). The power-law photon indices of spectra 1, 3, and 4 are tied, in 
accordance with Table [TJ 



